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1. Introduction 

The Fermilab Lattice and MILC collaboration's program includes calculations of the hadronic 
matrix elements for weak D and B meson decays, in particular, the decay constants /d, //)^, fs, 
and /b^ and the semileptonic form factors iov B ^ nlv , D ^ K{K)iv, and B ^ D*£v. In this 
work we present a perturbative matching calculation of the relevant current renormalizations to 
one-loop order. The numerical simulations for the above physics analyses use MILC ensembles 
with improved glue and 2+1 Asqtad staggered sea quarks The light valence quarks are also 
generated from Asqtad staggered quarks and converted to naive quarks. The heavy (charm and 
beauty) quarks are treated with the Fermilab action. See Ref. for more details on the actions 
and parameters used in the numerical simulations. 

2. Definitions 

In this work we follow the analysis of Ref. ||3|], where the one-loop corrections to heavy-light 
and heavy-heavy current renormalizations were calculated for Fermilab heavy and Clover light 
quarks with Wilson glue. 

The heavy-light currents have the form 

4"'''=Wk^^¥i, (2.1) 

where = 7^ or 7^75 and Yi denotes a naive Asqtad Dirac spinor. The Fermilab Dirac spinor, Yh, 
is rotated by 

Yh = ¥['^+adiY-^, (2.2) 

with the tree-level coefficient 

'^' = V^ ^ ^2-^^ 

2 + mofl: 2(1 +moa) 

The heavy-heavy currents have the form 

jf'^' = Y,r^Yh', (2.4) 

where now both spinors are rotated Dirac spinors. Since the heavy quarks are rotated, the lattice 
currents of Eqns. ( pH] ) and (2.4) include the leading order tree-level discretization corrections. 



The current renormalization is defined as 

,(1/2). 7(l/2)xcont 



^Jr- .7(1/2). 7(1/2).!,^ ' ^'■■^> 

where A/j. are the vertex corrections and Z2/, (Z2/) are the heavy (light) quark wave function renor- 
malizations. 

We factor out the dominant mass dependence due to the tree-level wave function renormaliza- 
tion of the heavy Fermilab quark by defining the perturbative expansion as 

.-'"'W^ = l+ggzJ!W + ..., (2.6) 
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where the heavy quark masses are defined as usual, 

mfc? = log(l +mofl') , moa = \ / {iKh) - \ / {iK^-ri,) . (2.7) 

Since Zy^ for degenerate masses is easy to calculate nonperturbatively, it is useful to define 

7/1/ 

In this case, the dominant mass dependence cancels by construction. 
Analogously, for heavy-heavy currents we have: 

^7(1/2)a 7(l/2)Ncont 
7/1/1' _ V^2/i ^^Jy^lh' ) /9 qx 

/7(l/2)., 7(l/2)yat ■ ^ ' 

Taking the leading mass dependence out again, the perturbative expansion is defined as 

.-('"^"'w'W^zf = 1+ggzf W + . (2.10) 

Finally, the p factors for heavy-heavy currents are defined as 

7/1/1' 

#-^%=^ = l+^^pf'^' + -- (2-11) 
'Z Z 



3. Procedure 



In this work we use the automated perturbation theory techniques developed by Liischer and 
Weisz ^ to generate the Feynman rules for the lattice actions. We then integrate the loop diagrams 
by "brute-force" using VEGAS [^]. The advantage of using automated perturbation theory is that 
it is relatively easy to switch actions [^]. Indeed, we have results for the current renormalizations 
for two gluon actions, two light quark actions and the heavy quark action. 

The one-loop diagrams for the vertex corrections (including the rotations) are given in Ref. [^]. 
We have performed the following tests of our calculation: 

• For the automated perturbation theory code, we have compared our vertices and propagators 
against known results. 

• We have written two independent programs for calculating the current renormalizations 
based on the automated perturbation theory code. 

• We have a third independent calculation of the current renormalizations using traditional 
semi-analytic methods. 

• Our results for the heavy-heavy currents agree with those of Ref. ^ when we switch from 
the improved gluon propagator to the Wilson gluon propagator. We also reproduce the results 
of Ref. for heavy-light currents with Clover light quarks and Wilson glue. 

• Our result for the Asqtad naive wave function renormalization agrees with Ref. [Q], and our 
result for the naive-naive vertex correction with Wilson glue agrees with Ref. 
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4. Results 



-0.02 



-0.04 



> 

N 




-0.06 - 



-0.08 



Figure 1: z'y^ for equal masses as a function of mf' . 

Figures 1^ illustrate our results for the Z's, p's, and g*'s as functions of the heavy-quark 
mass, nif^a. The <7*'s are calculated from the log moments using Eqn. (13) of Ref. [^. 

Our results for the heavy-heavy currents are very similar to those of Ref. [^], since they differ 
only in the gluon propagator. The main features of the mass dependence are the same. Figures 1-2 
show results for the degenerate mass V4 current. We also have results for the other currents (Vi, A4, 
A,) as well as results for currents with unequal masses. In the massless limit we find 

Z^'[^l(mfl =0) = -0.10056(3) , (4.1) 

in good agreement with Ref. [|Io|]. This is another test of our calculation. 

Figure 3 shows a comparison of the current renormalization of the heavy-naive A4 current with 
the corresponding p factor, and Figure 4 shows Py''^' and p^''^'. First, the general features of the 
heavy-quark mass dependence are similar to the results of Ref. [^]. Second, p^^'^' is significantly 
smaller than Z^^'^' over the relevant mass range. Hence, the cancellation between the numerator 



and denominator of Eq. (2.8) already observed in Ref. pp also takes place for heavy-naive currents. 
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Figure 2: q*a for Zy^ for equal masses as a function of ni\^ . 



In the massless limit we find: 



p^^W(mfl =0) = -3.038 (2) -10-3 



(4.2) 



p;:'''<(m';'' =0) = -3.05 (5) -10-3 
We also have results for the naive ^4 current renormalization. In the massless limit we find 

(4.3) 



Zy}^\mo = 0) = -0.10457 (4). 



We have studied the mass dependence of Zyj by varying mq between zero and the strange quark 
mass. We find that Zyj^^ is essentially independent of mo. 

In summary, we have calculated the current renormalizations relevant for the numerical analy- 
ses of heavy-light decay constants and semileptonic form factors performed by the Fermilab Lattice 
and MILC collaborations. We calculate the full mass dependence of the Z's and p's. The one-loop 
corrections to the p factors are small. They vary roughly between 0.4% and 4%, depending on 
lattice spacing. 
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Figure 3: Comparison of Z^^^ ' with p^^ 
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Figure 4: p factors for heavy-naive currents 
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